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Next-to-leading order corrections to jet cross sections in deep inelastic scattering at HERA are studied. The 
predicted jet rates allow for a precise determination of a 3 (/Xi?) at HERA over a wide range of fiR. We argue, that 
the "natural" renormalization and factorization scale is set by the average fcf? of the jets in the Breit frame and 
suggest to divide the data in corresponding < k§ > intervals. Some implications for the determination of the 
gluon density and the associated forward jet production in the low x regime at HERA are briefly discussed. 



1. Introduction 

Deep inelastic scattering (DIS) at HERA is a 
copious source of multi-jet events. Typical two- 
jet cross sections are in the 100 pb to few nb 
range and thus provide sufficiently high statistics 
for precision QCD tests [ Clearly, next-to- 
leading order (NLO) QCD corrections are manda- 
tory on the theoretical side for such tests. Full 
NLO corrections for one and two-jet production 
cross sections and distributions are now available 
and implemented in the ep — > n jets event gener- 
ator MEPJET, which allows to analyze arbitrary 
jet definition schemes and general cuts in terms 
of parton 4-momenta[ 0, §|. A variety of topics 
can be studied with these tools. 

1) The determination of cy s (/z.r) over a range 
of scales fin from dijet production at HERA: The 
dijet cross section is proportional to a s (fifi) at 
leading order (LO), thus suggesting a direct mea- 
surement of the strong coupling constant. How- 
ever, the LO calculation leaves the renormal- 
ization scale fiji undetermined. The NLO cor- 
rections substantially reduce the renormalization 
and factorization scale dependencies which are 
present in the LO calculations and thus reliable 
cross section predictions in terms of a s (mz) (for 
a given set of parton distributions) are made pos- 
sible. Clearly, a careful study of the choice of 
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scale in the dijet cross section is needed in or- 
der to extract a reliable value for a s (M z ). We 
will present some arguments and studies of the 
scale dependence of the NLO dijet cross section 
which suggest that the average of the jets in 
the Breit frame provides the "natural" scale for 
multi jet production in DIS. Here, (k^(jet)) 2 is 
defined by 2 E?(l — cos 0jp), where the subscripts 
j and P denote the jet and proton, respectively 
(all quantities are defined in the Breit frame). 

2) The measurement of the gluon density in 
the proton (via jg — > qq): The gluon density 
can only be indirectly constrained by an analy- 
sis of the structure function Fi at HERA [ 
The boson gluon fusion subproccss dominates the 
two jet cross section at low x and allows for a 
more direct measurement of the gluon density 
in this regime. A first LO experimental anal- 
ysis has been presented in [ ||. NLO correc- 
tions reduce the factorization scale dependence 
in the LO calculation (due to the initital state 
collinear factorization, which introduces a mix- 
ture of the quark and gluon densities according 
to the DGLAP (Altarelli-Parisi) evolution) and 
thus reliable cross section predictions in terms 
of the scale dependent parton distributions are 
made possible. Some implications for the deter- 
mination of the gluon density have been discussed 
in[|. 

3) Associated forward jet production in the 
low x regime as a signal of BFKL dynamics: 
BFKL evolution [ ^| leads to a larger cross sec- 
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tion for events with a measured forward jet (in 
the proton direction) with transverse momentum 
p l f b (j) close to Q than the DGLAP [ evolu- 
tion. Clearly, next-to-leading order QCD correc- 
tions for fixed order QCD, with Altarelli-Parisi 
(DGLAP) evolution, are mandatory on the theo- 
retical side in order to establish a signal for BFKL 
evolution in the data. We discuss these correc- 
tions in section 3. 

The importance of higher order corrections and 
recombination scheme dependencies of the two 
jet cross sections for four different jet algorithms 
(cone, kr [ ||, JADE, W) was already discussed 
in [ ||, || . While the higher order corrections and 
recombination scheme dependencies in the cone 
and kx schemes are small, very large corrections 
appear in the W^-scheme. We conclude from these 
studies that the cone and fcy schemes appear bet- 
ter suited for precision QCD tests and will con- 
centrate only on those schemes in the following. 

The goal of a versatile NLO calculation is to 
allow for an easy implementation of an arbitrary 
jet algorithm and to impose any kinematical res- 
olution and acceptance cuts on the final state 
particles. This is best achieved by performing 
all hard phase space integrals numerically, with 
a Monte Carlo integration technique. This ap- 
proach also allows an investigation of the recom- 
bination scheme dependence of the NLO jet cross 
sections. For dijet production at HERA such a 
NLO Monte Carlo program is MEPJET[ §. The 
calculation is based on the phase space slicing 
method and on the technique of universal "cross- 
ing functions" [ ^, |l0|, jll). More details are de- 
scribed in [ |^] and we do not repeat them here. 
This technique can also be extended to the case 
of massive quark production [ |l2| . An alternative 
technique for the calculation of NLO corrections 
in jet physics is the "subtraction" method [ |l3|. 

2. Choice of the Renormalization and Fac- 
torization Scale in Multijet Production 

Jet production in DIS is a multi-scale problem 
and it is not a priori clear at which scale a s is 
probed. In the following we will discuss the choice 
of the renormalization and factorization scale for 
dijet production at HERA. As mentioned before, 



the NLO corrections are expected to reduce the 
scale dependencies in the LO predictions provided 
the scale is of the order of the typical hardness of 
the partonic process. 

The following studies are done for the cone 
algorithm (which is defined in the laboratory 
frame) and the distance AR = sj (Arf) 2 + (A<^>) 2 
between two partons decides whether they should 
be recombined to a single jet. Here the variables 
are the pseudo-rapidity r/ and the azimuthal an- 
gle <fi. We recombine partons with AR < 1. Fur- 
thermore, a cut on the jet transverse momenta 
of pt{j) > 5 GeV in the lab is imposed. We 
employ the two loop formula for the strong cou- 
pling constant with a value for A^- consistent 
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with the value from the parton distribution func- 
tions. The value of a s is matched at the thresh- 
olds fin — m q and the number of flavors is fixed 
to nf = 5 throughout, i.e. gluons are allowed 
to split into five flavors of massless quarks. A 
running QED fine structure constant a{Q 2 ) is 
used. In addition the following set of kinemat- 
ical cuts is imposed on the initial virtual photon 
and on the final state electron and jets: We re- 
quire 40 GeV 2 < Q 2 < 2500 GeV 2 , 0.04 < y < 1, 
an energy cut of E(e') > 10 GeV on the scat- 
tered electron, and a cut on the pseudo-rapidity 
Tj = — lntan(6'/2) of the scattered lepton and jets 
of \<q\ < 3.5. 

First studies of the scale dependence of the di- 
jet cross section in the cone scheme are presented 
in [ E|. We have considered scales related to the 
scalar sum of the parton transverse momenta in 
the Breit frame, JV PtM, an( i the virtuality Q 2 
of the incident photon. In the following we will 
also consider scales related to JV k^(i). Here, 
(k^(i)) 2 is defined by 2E 2 (1 — cos^p), where 
the subscripts i and P denote the final parton (or 
jet) and proton, respectively (all quantities are 
determined in the Breit frame). The Breit frame 
is characterized by the vanishing energy compo- 
nent of the momentum of the exchanged photon. 
Both the photon momentum 

q= (0,0,0, -2xE), -q 2 = Q 2 = Ax 2 E 2 (1) 

and the proton momentum 

P = E(1, 0,0,1) (2) 
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are chosen along the z-direction. x is the stan- 
dard Bjorken scaling variable. In the parton 
model, the incoming quark with momentum p = 
xE(l, 0, 0, 1) collides clastically with the virtual 
boson and is scattered in the opposite direc- 
tion (the 'current hemisphere') with momentum 
p' = xE(l, 0,0,-1); therefore, (fcf(p')) 2 = Q 2 
in the limit of the quark parton model, whereas 
pf(p') = 0. 

The kinematics for dijet production is more 
complex: the momentum fraction rj of the incom- 
ing parton must be larger than x since rrijj = 
Q 2 (r]/x — 1) and, in general, the jets have a non- 
vanishing transverse momentum with respect to 
the 7*-proton direction. At LO, i.e. for mass- 
less jets, the relation between k^U) and p^U) = 



Ej sin ( 



eads: 



jP r ' 



PtU) 




cos 9 



(3) 



Obviously, k^U) > p^, and one can also show 
that J2j ^t(j) > Q- Thus, J^j ^t(j) is approxi- 
mately given by the harder of the two scales Q 
and J^j PtU) I @- ^or l ar S e dijet invariant 
masses (i.e. for "true" two jet kinematics) one has 
i) » x, the dijet system will be strongly boosted 
in the proton direction and k^U) w PtU) >> 
Q/2, as long as very small scattering angles in the 
center of mass frame are avoided. For rrijj ~ Q 
or smaller, on the other hand, J^j ^t(j) ~ Q 
and, typically, both are considerably larger than 
PtU) < 171 jj ( m particular for Q >> rrijj, 
which corresponds to the parton model limit). 

Thus J^j ^t GO smoothly interpolates between 
the correct limiting scale choices, it approaches 
Q in the parton limit and it corresponds to the 
jet transverse momentum when the photon vir- 
tually becomes negligible. It appears to be the 
"natural" scale for multi jet production in DIS. 

In [ D we found that the scale dependence of 
the dijet cross section does not markedly improve 
in NLO for /i 2 = £,Q 2 . This is also shown in Fig. 1 
(dotted curves) where the dependence of the two- 
jet cross section on the scale factor £ is shown. We 
used the parton distribution functions set MRS 
D-' [ [15) and employed the two loop formula for 
the strong coupling constant both in the LO and 
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Figure 1. Dependence of the two-jet exclusive 
cross section in the cone scheme on the scale fac- 
tor £. The dashed curves are for n 2 R = /i 2 ^ = 
f (Ei Pt (0) 2 - Choosing (£ { fe| (i)) 2 as the ba- 
sic scale yields the solid curves. Choosing Q 2 as 
the basic scale yields the dotted curves. Results 
are shown for the LO (lower curves) and NLO 
calculations. 



NLO curves. For scales related to J^i Pt(0 the 
uncertainty from the variation of the renormaliza- 
tion and factorization scale is markedly reduced 
compared to the LO predictions (dashed curves 
in Fig. 1). Here £ is defined via 



(J? F 



(4) 



The resulting £ dependence for fi R 
£ E» (*)) 2 1S shown as the solid lines in Fig. 1. 
In this case, the NLO two-jet cross section is es- 
sentially independent on £ for £ < 2. Hence, the 
theoretical uncertainties due to the scale variation 
are very small suggesting a precise determination 
of a s (< k§ >) for different < kj, > bins, where 



J=l,2 



(5) 



Fig. 2 shows the < fc|? > distribution for the 
NLO exclusive dijet cross section. We used the 
parton distribution functions set GRV [ 16 and 
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Figure 2. 

jet exclusive cross section. The parameters are 
explained in the text. 



< kj, > distribution for the two- 
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Figure 3. < pf > distribution for the three 
< k^ > bins shown in Fig. 1. 
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Figure 4. Q 2 distribution for the three < k^ > 
bins shown in Fig. 1. 



Table 1 

Jet cross sections in pb for the three < > bins. 



bin 1 bin 2 bin 3 



GRV£=l/4 497 

GRV£=1/16 504 

GRV £ = 1 488 

MRSD-' £ = 1/4 487 



320 
306 
322 
322 



146 
134 
151 
148 



In addition to the 

B 



& = Ml = 1/4 (E 

cuts imposed in Fig. 1, we require p^, > 4 GeV 
for each jet and Ej C?) > 10 GeV. We have di- 
vided the NLO cross section in the following three 



< k^ > bins. 



bin 1 
bin 2 
bin 3 



5 GeV < < jfef > < 10 GeV 
10 GeV < < fc| > < 15 GeV 
15 GeV < < fc| >. 
Table 1 shows the corresponding NLO cross sec- 
tions for these three bins for two different sets of 
parton distributions and for different values for 
the scale factor £. The theoretical uncertainties 
for the NLO dijet cross section from these num- 
bers are very small, in particular for the first two 
bins. Fig. 3 and 4 show the < pj, > and Q 2 distri- 
bution for the NLO exclusive dijet cross section 
for these three bins. Whereas the < pj, > and 
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Q 2 distributions are fairly similar to the < fcf > 
distribution for the lowest bin, large differences 
are found for the other bins. This partly explains 
the rather different scale dependence observed in 
Fig.l. 



3. Forward Jet Production in the Low x 
Regime 

Deep inelastic scattering with a measured for- 
ward jet with relatively large momentum fraction 
Xj e t (in the proton direction) and p 2 } ab (j) ~ Q 2 
is expected to provide sensitive information about 
the BFKL dynamics at low x [ [l7|, |l8| . In this re- 
gion there is not much phase space for DGLAP 
evolution with transverse momentum ordering, 
whereas large effects are expected for BFKL evo- 
lution in x. In particular, BFKL evolution is ex- 
pected to substantially enhance cross sections in 
the region x « Xj et [ p|. In order to ex- 
tract information on the ln(l/a;) BFKL evolution, 
one needs to show that cross section results based 
on fixed order QCD with DGLAP evolution are 
not sufficient to describe the data. Clearly, next- 
to-leading order QCD corrections to the DGLAP 
predictions are needed to make this comparison 
between experiment and theory. 

In Table | we show numerical results for the 
multi jet cross sections with (or without) a for- 
ward jet. The LO (NLO) results are based on the 
LO (NLO) parton distributions from GRV [ [jU 
together with the one-loop (two-loop) formula for 
the strong coupling constant. Kinematical cuts 
are imposed to closely model the HI event se- 
lection [ |l^]. More specifically, we require Q 2 > 
8 GeV 2 , x < 0.004, 0.1 < y < 1, an energy cut 
of E(e') > 11 GeV on the scattered electron, and 
a cut on the pseudo-rapidity rj — — lntan(#/2) of 
the scattered lepton of -2.868 < rj(e') < -1.735 
(corresponding to 160° < 6(1') < 173.5°). Jets 
are defined in the cone scheme (in the laboratory 
frame) with AR = 1 and \r)(j)\ < 3.5. We re- 
quire a forward jet with Xj e t = p z (j)/Ep > 0.05, 
E(j) > 25 GeV, 0.5 < p 2 T {j)/Q 2 < 4, and a cut 
on the pseudo-rapidity of 1.735 < r/(j) < 2.9 (cor- 
responding to 6.3° < 0(j) < 20°). In addition all 
jets must have transverse momenta of at least 4 
GeV in the lab frame and 2 GeV in the Breit 



Table 2 

Cross sections for n-jet exclusive events in DIS at 
HERA. See text for details. 



with 
forward jet 



without 
forward jet 



1 jet (LO) 

2 jet (LO) 

2 jet (NLO) 

3 jet (LO) 



pb 
19.3 pb 
68 pb 
30.1 pb 



9026 pb 
2219 pb 
2604 pb 
450 pb 



frame. 

The cross sections of Table || demonstrate first 
of all that the requirement of a forward jet with 
large longitudinal momentum fraction (xj e t > 
0.05) and restricted transverse momentum (0.5 < 
PtU)/Q 2 < 4) severely restricts the available 
phase space, in particular for low jet multiplic- 
ities. The 1-jet exclusive cross section vanishes 
at LO, due to the contradicting x < 0.004 and 
Xj e t > 0.05 requirements. For x << Xj e t, a high 
invariant mass hadronic system must be produced 
by the photon-parton collision and this condition 
translates into 



2E(j)rriT e 



£7, parton 
\ X 



l »Q 2 



(0) 



where ttit and y are the transverse mass and ra- 
pidity of the partonic recoil system, respectively. 
Thus a recoil system with substantial transverse 
momentum and/or invariant mass must be pro- 
duced and this condition favors recoil systems 
composed out of at least two additional energetic 
partons. 

As a result one finds very large fixed order per- 
turbative QCD corrections (compare 2 jet LO and 
NLO results with a forward jet in Table ||). In 
addition, the LO (0(a 2 )) 3-jet cross section is 
larger than the LO (0(a s )) 2-jet cross section. 
Thus, the forward jet cross sections in Table |^ 
are dominated by the (0(a 2 )) matrix elements. 
The effects of BFKL evolution must be seen and 
isolated on top of these fixed order QCD effects. 
We will analyze these effects in a subsequent pub- 
lication. 
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4. Conclusions 

The calculation of NLO perturbative QCD cor- 
rections has received an enormous boost with the 
advent of full NLO Monte Carlo programs [ ||, ||o| . 
For dijet production at HERA the NLO Monte 
Carlo program MEPJET [ |j allows to study jet 
cross sections for arbitrary jet algorithms. Inter- 
nal jet structure, parton/hadron recombination 
effects, and the effects of arbitrary acceptance 
cuts can now be simulated at the full 0(a 2 s ) level. 
We found large NLO effects for some jet defini- 
tion schemes (in particular the VF-scheme) and 
cone and kx schemes appear better suited for pre- 
cision QCD tests. 

The extraction of gluon distribution functions 
is now supported by a fully versatile NLO pro- 
gram. Preliminary studies show that large NLO 
corrections are present in the Bjorken x distribu- 
tion for dijet events, while these effects are mit- 
igated in the reconstructed Feynman x (xi) dis- 
tribution, thus aiding the reliable extraction of 
9(xi,lJ, 2 F ). 

For the study of BFKL evolution by consid- 
ering events with a forward "Mueller" -jet very 
large QCD corrections are found at 0(a 2 s ). These 
fixed order effects form an important background 
to the observation of BFKL evolution at HERA. 
They can now be studied systematically and for 
arbitrary jet algorithms. 
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